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We present an explicit model where the decay of an _R-parity even scalar S with 0{TeV) mass 
is the origin of non-thermal dark matter. The correct relic abundance can be produced for both 
large and small annihilation rates in accordance with the Fermi constraints on the annihilation 
cross-section. This scenario has advantages over that of non-thermal dark matter from modulus 
decay. First, branching ratio for production of i?-parity odd particles can be made quite small by 
a combination of 5* couplings to matter fields and kinematic suppression, enabling us to obtain the 
observed dark matter relic density in cases of thermal underproduction as well as overproduction. 
Second, gravitino production is naturally suppressed by the virtue of decaying scalar belonging to 
the visible sector. The decaying scalar can also successfully generate baryon asymmetry of the 
universe, and may provide an explanation for the baryon-dark matter coincidence puzzle. 



I. INTRODUCTION 

Weakly interacting massive particles (WIMPs) are 
promising dark matter (DM) candidates that can explain 
the DM relic abundance, as precisely measured by cos- 
mic microwave background experiments pi, via thermal 
freeze-out of annihilation in the early universe. The nom- 
inal DM annihilation rate for this scenario, called the 
WIMP miracle, is {as.nnv) « 3 x 10"^^ cm^ s"^. 

WIMPs typically arise in models of particle physics 
beyond the standard model (SM). In supersymmetric 
(SUSY) models with conserved i?-parity, the lightest su- 
persymmetric particle (LSP) is the DM candidate. The 
lightest neutralino is the most suitable candidate with 
prospects for detection in various direct and indirect 
searches. However, in large regions of SUSY parameter 
space the thermal relic abundance of the LSP is differ- 
ent from the observed value. Typically, if the lightest 
neutralino is a Higgsino or Wino, the annihilation cross- 
section is large compared to the nominal value, while for 
a Bino the annihilation cross-section is small. 

A mainly Higgsino LSP is motivated by considerations 
of naturalness [2HZj ■ If the Higgsino mass is in the sub- 
TeV region, the annihilation rate is typically larger than 
3 X 10~^^ cm^ s~^, thus resulting in an insufficient ther- 
mal relic abundance. The Fermi-LAT data ^ constrains 
the DM annihilation rate at the present time. These 
bounds do not allow much room for too large an annihi- 
lation cross-section. Further, for lower values of the LSP 
mass, Fermi-LAT data appears to prefer smaller annihi- 
lation cross-sections. This prefers Bino type LSP, which 
will have a thermally overproduced abundance. 

Therefore, since the presently motivated DM annihila- 
tion cross-sections are either larger or smaller than the 
nominal thermal freeze-out value, one is naturally led to 
consider non-thermal scenarios to obtain the correct DM 
relic density. A late-decaying field S that reheats the 
universe below the freeze-out temperature Tf ^ 



{x denoting the DM particle) can provide a non-thermal 
origin for DM. If the annihilation rate is larger than the 
nominal value, the LSPs produced from S decay undergo 
residual annihilation before reaching their final abun- 
dance. Models that prefer Higgsino DM can be accom- 
modated in such a scenario. If the number density of 
LSPs produced from S decay is very small and/or the 
annihilation rate is smaller than the nominal value, an- 
nihilation will be inefficient and the final DM relic density 
will the same as that from S decay. This scenario can ac- 
commodate both types of models with Higgsino and Bino 
DM. 

A gravitationally coupled modulus field, typically aris- 
ing in SUSY and superstring-inspired models [9 , serves 
as a standard candidate for the field S. Moduli that are 
heavier than 50 TeV decay before the onset of big bang 
nucleosynthesis (BBN), and hence will not ruin its suc- 
cessful predictions on the primordial abundance of light 
elements. The decay also provides a non-thermal origin 
for DM production [TUIIII]. However, in such a scenario, 
the branching ratio for gravitino production from late de- 
cay (denoted by Br3/2) must be sufficiently suppressed. 
Otherwise, the decay of gravitinos will lead to DM over- 
production or, for gravitino mass below 40 TeV, will ruin 
the success of BBN. It is not easy to satisfy this require- 
ment in scenarios where modulus decay is the source of 
DM production. For example, in simple models such as 
KKLT [l^ Br3/2 turns out to be too large by a factor of 
~ 10^ HI]. Furthermore, it is difficult to obtain a sufh- 
ciently small DM abundance from moduli decay without 
residual annihilation. Even if one suppresses LSP pro- 
duction from two-body decays of moduli, the three-body 
decays will still be unacceptably large |13j . Suppressing 
the total decay rate of the modulus can help the situa- 
tion, but it involves non-trivial conditions on the Kahler 
geometry of the underlying effective supergravity theory. 

In this paper we explore the late decay of a visible sec- 
tor scalar field S that serves as the origin of non-thermal 
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DM. It has the foUowing advantages: (1) Gravitino pro- 
duction from S is naturally suppressed for gravitino mass 
as low as 0{TeV) by virtue of the fact that S" is a vis- 
ible sector field and (2) Due to kinematic suppression, 
i?-parity odd particles are produced at one-loop or two- 
loop level, which can yield Br^^ <^ 10~^. This immedi- 
ately opens up the option of obtaining non-thermal DM 
purely from branchings of the scalar field, without un- 
dergoing further annihilation. As we discussed, for Bino 
LSP, this is the only option. For Higgsinos, this allows 
both options of a relic density set by pure branching ra- 
tios, as well as by annihilation. 

We present an explicit model where a visible sector 
scalar S with a mass ms ~ 0{TeV) and i?-parity charge 
-1-1 has direct couplings to new colored fields X, X . If S* 
is lighter than X, X particles, as well as all colored su- 
perparticles, it dominantly decays into two gluons at one- 
loop level. S decay also produces DM particles via loops, 
but this mode is suppressed by powers of m^/ms and 
51,2/53, with gi^2,3 being the C/(l)r, SU{2)w, SU{3)c 
gauge couplings respectively. As we will see, this model 
can yield the correct DM relic abundance for of the ther- 
mal underproduction and overproduction cases, and also 
lead to successful late-time baryogenesis. 

The paper is organized as follows. In Section II, we 
briefly review non-thermal DM form modulus decay and 
the associated problems. In Section III, we introduce 
a model with late decaying scalar field and discuss its 
properties. In Section IV, we show that this model can 
successfully explain the DM content for both large and 
small annihilation cross sections. In Section V, we dis- 
cuss baryogenesis in this model and point out that the 
presented non-thermal DM scenario can also address the 
baryon-DM coincidence puzzle. We conclude the paper 
in Section VI. 



we find the necessary condition for S to be dominant at 
the time of decay 
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Decay of S reheats the universe to a temperature Tj. ^ 
(FsAfp)^/^. As a numerical example, for ms C(TeV) 
and Ti- ~ 3 MeV (in order to be compatible with BBN), 
Eq. Q imphes S dominance for Sq » 10^^ GeV. 

If S dominates the universe at a temperature Tdom, we 
will have 



Pr, after 
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(2) 



where "before" and "after" are in reference to the epoch 
of S decay, and we have used the fact that Pr^aftcr — Ps- 
It is seen from Eq. ^ that S decay releases a large 
entropy that dilutes any pre-existing quantity in the 
thermal bath. For the above numerical example where 
ms ~ ©(TeV) and Tr ~ 3 MeV, the entropy release fac- 
tor can be as large as 10®. 



B. Dark matter from late decay 



Provided that < Tr 



^/25, decay of S will di- 



lute any thermally produced DM by a large factor as 
mentioned above. However, S decay itself produces DM 
particles. The abundance of non-thermally produced DM 
is given by 



II. NON-THERMAL DARK MATTER FROM A 
LATE DECAYING SCALAR FIELD 

In this section, we discuss the various possibilities for 
non-thermal production of DM from a late decay. 



A. Reheating by late decay 

We consider a scalar field S with mass ms and de- 
cay width F5. Assuming that S has acquired a large 
vacuum expectation value during inflation, it will start 
oscillating about the minimum of its potential with an 
initial amplitude when the Hubble expansion rate is 
H ~ ms- Oscillations of S behave like matter, with an 
initial energy density ps — to|.S'o/2. The energy density 
of the universe at this time, dominated by thermal bath, 
is = 3m|A/p. 

The quantity psi P-c is rcdshifted cx a, with a being the 
scale factor of the universe. After using the fact that iJ 
is redshifted oc a^^ for a radiation-dominated universe. 



Ys Br 



thr 



(3) 



Here I5 = ZT^jAms^ Br^ is the branching fraction for 
production i?-parity odd particles from S decay, and 
(n^/s)thr) denotes DM abundance obtained via thermal 
freeze-out that is related to the observed DM relic abun- 
dance (n^/s)obs through: 



(^) 



thr 



3 X 10~2S cm^ s-i 



S / obs 



V S / obs 



5 X 10 



-10 



(CannW) 



1 GeV\ 



(4) 
(5) 



The abundance of DM particles immediately after 
their production from S decay is given by YgBr^. If 
n^^Os.miV) < H{T-c), DM annihilation will be inefficient 
at temperature In this case, the final DM relic abun- 
dance will be given by the first inside the brackets in 
Eq. (j3|. On the other hand, if n^{(j^^^v) > H{Tr), an- 
nihilation will be efficient right after S decay. This will 
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somewhat reduce the abundance of DM particles pro- 
duced from S decay, in which case the final relic density 
will be given by the second term inside the brackets in 
Eq. (|3|. 

There are therefore two possible scenarios for obtaining 
the correct DM relic density from S decay: 

• Annihilation Scenario: If (aann) > 3 x 
10^^^ cm'^ s^^, then (n^/s)thr < {nchi/s)ohs (hence 
"thermal underproduction"). The large annihila- 
tion cross section can reduce the abundance of DM 
particles produced from S decay to an acceptable 
level, provided that: 



T, = Tt 



3 X 10^26 pj^^3 g-1 



The final DM abundance will then be given: 



(-) 

V S /thr 



3 X 10^26 g-1 / jj^ 



(O'annl') 



(6) 



(7) 



This scenario can work well in the case of Higgsino 
DM, for which (ffannw) > 3 X IQ-^^ cm^ s"^) as 
mentioned before, provided that the reheat tem- 
perature from S decay satisfies Eq. 

• Branching Scenario: If Eq. (|6| is not satisfied, 
then annihilation will be rendered ineffective. This 
happens if is too low and/or ((Tannw) is too small. 

The first possibility is that (cann > 3 x 
10~^^ cm'' s""'^, but Tr is lower than that given 
in Eq. In this case non-thermal Higgsino DM 
must be produced via "Branching Scenario" . 

On the other hand, we note that Eq. (l6| can never 
be satisfied if (o-annw) < 3 x 10~^^ ciir s~^. It is 
seen from Eq. Q that this results in (n^/s)thr > 
{n-^/s)ohs (hence "thermal overproduction"). This 
leaves "Branching Scenario" as the only possibility 
for non-thermal DM production in this case. Bino 
DM provides a prime example of this case. 

The final DM abundance will be the same as that 
produced from S decay, which follows 



"Annihilation Scenario" can be realized in this frame- 
work. We have recently discussed non-thermal Hig- 
gsino DM from modulus decay in scenarios with mixed 
anomaly-modulus mediation of supersymmetry break- 
ing 5 . For a typical modulus mass ms ^ 0(1000) 
TeV'in this scenario, one finds T, ~ ©(GeV). Eq. ^ 
can be satisfied for Higgsino annihilation cross sections 
that are compatible with Fermi bounds [5] with a mass 
~ 100- 1000 GeV 0. 

"Branching Scenario" , however, is not easily realizable. 
The fact that Ys > 10~'^ requires Br^ < 5 x 10-*^ for 
> 100 GeV. Such a small Br^ may be obtained for 
two-body decays of the modulus [TD] , but three-body de- 
cays will inevitably set a lower bound Br-^ > 10"'^. Ob- 
taining the correct DM relic abundance within "Branch- 
ing Scenario" then requires that Ys be further lowered, 
which may happen by means of geometric suppression 

An additional challenge is posed by gravitino pro- 
duction from modulus decay, whose abundance follows 
(713/2/5) = Ys Br3/2, where Br3/2 is the branching frac- 
tion for gravitino production from modulus decay. Grav- 
itinos decay much later than the modulus. For exam- 
ple, in the scenario discussed in [5] m3/2 ^ 0(100) TeV, 
which implies they decay just before the onset of BBN. 
Gravitinos produce DM particles upon decay, which re- 
quires that 



< 5 X 10-1" 



1 GeV 

JTiv 



(9) 



For > 100 GeV, this results in the bound 

("3/2/5) < 5 X IQ-i^. Since Ys ^ lO^*", we then need 
Br3/2 < 10"^. However, in the simples example based 
on KKLT model [12] we have Br3/2 10^^. One may 
lower Br3/2 by modifying the Kahler potential, but a suc- 
cessful implementation that does not affect other aspects 
of the non-thermal DM scenario is a non-trivial task. 

To summarize, a completely successful non-thermal 
DM scenario from modulus decay is challenging, in par- 
ticular in models with thermal overproduction (notably 
Bino DM). The reason being that Br^ and Br3/2 are typ- 
ically too large in this case. 



"a Brv. 
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III. LATE DECAY OF A VISIBLE SECTOR 
SCALAR 



C. Challenges for non-thermal dark matter from 
modulus decay 

Modulus decay provides a natural scenario for non- 
thermal DM jlOj. Moduli heavier than 50 TeV decay 
before the onset of BBN, which allows their utilization 
as the source of DM production. The modulus decay rate 
is given by Ts ~ m|/27rM| [B]. For ms > 50 TeV, this 
results in Ys > 10^'^. 



In this section, we present an explicit model of non- 
thermal DM from late decay of a visible sector field. We 
describe the model and its field content, and show how 
it can address the abovementioned issue with Br3/2- 

A. The Model 

The visible sector consists of the minimal supersym- 
metric standard model (MSSM) augmented with extra 
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superfields: 

(i) A singlet S whose decay is the origin of DM. 

(m) Two flavors of iso-singlet color triplets ^1,2, -?i,2 
with hypercharges +4/3, —4/3 respectively. 

(in) Two flavors of singlets A^i^2- 

Multiple flavors of X, X, N are introduced to accom- 
modate late-time baryogenesis (which we will discuss in 
Section V). 

We choose charge assignments under i?— parity such 
that i?-parity is conserved, hence LSP is stable. This 
implies that X, X, S and N are i?-parity even, while 
their SUSY partners are i?-parity odd. We also assume 
that the lightest MSSM neutralino is the DM candidate. 

The superpotential of the visible sector is Wvisibic = 
W'mssm + Wn,x + Ws, where 



Wn.x ^ X'd^d^X + XNu^X 



NN + mxXX , 



and 



Ws = hSXX + ]-MsS 



(10) 



(11) 



For simplicity, we have omitted the flavor and color in- 
dices. Henceforth, we use the same symbol for superfields 
and their corresponding i?-parity even component fields, 
while i?-parity odd fields are distinguished by a". 

We note that some terms that are gauge-invariant are 
absent in Ws, namely SNu", N^, SNN, S^. The 
first two terms are forbidden by i?-parity, while invok- 
ing some other discrete or continuous symmetry (like an 
i?-symmetry) may help forbid or suppress the last two 
terms. 

Soft SUSY breaking terms typically make scalar com- 
ponents of chiral superfields heavier than their fermionic 
counterparts. In the case of S superfield, one has m| = 
M| -I- m^Qfj ± BMs, where mgoft denotes the soft mass of 
S and B is the _B-term associated with the superpoten- 
tial mass term for superfield S. We assume the following 
mass condition 



iriN < ms < mx 
ms > 2m^ . 



(12) 



It is also reasonable to assume that all i?-parity odd col- 
ored particles have a mass larger than ms- The rea- 
son being that soft breaking masses of these particles are 
driven toward large values at low energies by radiative 
corrections with SU{3)c gauge interactions while 5 is a 
singlet. 

These conditions imply that S cannot decay to either 
of X, X OT X, X fields. Moreover, it cannot decay to any 
_R-parity odd colored fields. Similarly, decay to i?-parity 



odd scalars is kinematically blocked since these masses 
are governed by the soft terms, which supposedly have 
the smallest value for a singlet field. S can however decay 
to DM particles, which is an essential part of this model. 

As a consequence, S dominantly decays into two gluons 
through the one-loop diagram shown in Figure [l]i The 
corresponding decay is given by 



S^gg 



2 X 



0.17 

"8^ 



47r2 



/ ms 



ms 



(13) 



Here 33 ^ 1 is the SU{3)c gauge couphng constant and 
the two flavors of X are taken into account. The precise 
expression for Ts^gg is given in the Appendix. 

The fact that S decay is loop suppressed combined 
with ms ^ mx help us obtain a sufficiently low reheat 
temperature Tj. (rsMp)^/^ < Tf. For example, for 
ms ^ 1 TeV and mx ^ 10 - 100 TeV, we find T, - 
©(GeV) if /i~ 10-^- 10-6. 

The magnitude of h needed is similar to the electron 
Yukawa coupling, and also comparable to typical values 
of neutrino Dirac Yukawa couplings in TeV scale see- 
saw models. Its smallness may be explained in different 
ways. There may exist some symmetry, broken at the 
scale of grand unified theories (GUT), under which S is 
charged. Then the term SXX arises from a higher-order 
Mp suppressed operator after symmetry breaking, and 
its strength will be suppressed by powers of Mqut/Mp. 
Moreover, see Eq. (13), we notice that the combination 
h/mx appears in F^, which acts like an effective cou- 
pling. Therefore one can make h larger by simultane- 
ously increasing mx with the same factor. For example, 
we find h - 0(1) if mx ^ 10" GeV. 

We also note that for a gravitationally interacting field 



S the decay rate will be oc m^/Mp. It is seen from 
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Eq. ([13| that S decay to gluons occurs with a strength 
much larger than that for a gravitational decay as long 
as AiT^mx/h <C Mp. This is clearly the case for 777.5 ~ 
O(TeV) and > 3 MeV, which justifies S belonging 
to the visible sector rather than a hidden sector with 
gravitationally suppressed coupling to matter. 




^ S can also decay into four-body final states containing 
(TV, u'^, d'^) through mediation of ofl^shell X, X. However, these 
decays are suppressed compared to the two-body decays by phase 
space factors and additional powers of ms/mx, and hence can 
be neglected. There are also subdominant decay modes of S that 
proceed via loop diagrams and are important for DM production 
and baryogenesis. We will discuss these modes later. 
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FIG. 1: Diagram showing the dominant one-loop decay of S. 
There is a similar diagram with X , X running in the loop. 



B. Gravitino Production 

Gravitinos can be produced from S decay if 7713/2 < 
nis- Considering that ms ~ 1 TeV, this is kineniatically 
possible only if 7713/2 <, C(TeV. This implies that there 
will be no gravitino production form S decay for a very 
wide mass range 7113/2 > ©(TeV). 

When kinematically possible, the important decay 
mode is S ^ G + S, whose decay width is given by [2] 



production of Higgsino and Bino DM from visible sec- 
tor decay. Recently, Higgsino type LSP has attracted 
significant attention from the natural SUSY perspective 
[2H7]. However, if the Higgsino mass is in the sub- TeV 
region, the annihilation rate is larger than the nominal 
value 3 x 10~^^ cm'^ s""'^, which yields insufficient thermal 
relic abundance. The concern with this scenario is that 
the current constraint from the Fermi-LAT data does not 
allow much room for the annihilation cross-section to be 
too large. These results appear to prefer smaller anni- 
hilation cross-sections for smaller values of LSP mass. 
This prefers Bino type LSP, which will have a thermally 
overproduced abundance. 



r 



S^G+S ^,j2- 



(14) 



Gravitinos with ©(TeV) mass decay long after BBN. 
Their abundance should be low enough in order not to 
ruin successful predictions of BBN. Since TO3/2 < 1715, 
and S is assumed to be lighter than all colored super- 
particles, gravitino decay modes can only be radiative. 
Successful BBN in this case requires that [T5] 



%/2 



< 10" 



-12 



(15) 



This also insures that gravitino decay will not overpro- 
duce DM particles with a mass ^ 500 GeV. 

The abundance of gravitinos produced from S decay 
follows 



"3/2 
s 
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(16) 



Branching Fraction for Decay to Dark Matter 
Particles 



The decaying scalar S has i?-parity charge +1, which 
implies that it can only decay to an even number of R- 
parity odd particles. If it was i?-parity odd instead, all 
of the decay modes would produce i?-parity odd parti- 
cles, thus resulting in Br^ 1. However, one can now 
obtain a small Br,^ if S decay to i?-parity odd particles 
is suppressed relative to that of S decay to gluons. 

The dominant mode for producing i?-parity odd par- 
ticles is 5" — >■ BB that proceeds through the one-loop 
diagram shown in Figure [2] (top diagram) . This mode 
is kinematically allowed if Bino is the DM particle, see 
Eq. ( 12 ). The corresponding decay width is given by (for 



more details see the Appendix) 



2 X 



0.12 

"8^ 



167r2 



^ ' ms,(17) 
\mx. 



For 7115 ~ O(TeV) the hmit in Eq. ^ is satisfied pro- ^j^^^^ ^(^^^ ^^^^^ coupling constant. Both 



vided that Tj. > 2.5 MeV. We see that gravitino produc 
tion will not be a problem for the entire range of Tj. al- 
lowed by BBN (i.e., > 3 MeV). Moreover, this condi- 
tion will be irrelevant altogether if 1173/2 > 1115, in which 
case S decay to the gravitino will be forbidden kinemat- 
ically. 

The fact that that gravitino overproduction from S de- 
cay is avoided comes as a direct consequence of S mainly 
decaying into visible sector fields (i.e., gluons), while its 
decay to gravitinos is Mp suppressed. The situation is 
very different when S* is a modulus field, as pointed out 
earlier, since all decay modes are gravitationally sup- 
pressed in that case. 



IV. DARK MATTER FROM VISIBLE SECTOR 
DECAY 

In this section we show how the model presented in the 
previous section can produce the observed DM relic abun- 
dance via the "Branching Scenario" and "Annihilation 
scenario" both. We also discuss successful non-thermal 



flavors of X have been taken. After using Eq. (13 1, this 
results in the following branching fraction 



Br, 



°- (Mr) 



(Bino LSP). (18) 



Inserting numerical values for the parameters 7115 '-^ 1 
TeV, - 100 GeV, 51 - 0.3, 53 ~ 1, and Yx = 4/3, 
we find Br^ ^ 1.1 x 10^^. This is much smaller than 
what one finds from modulus decay Br^^, > 10~^ [13] . 

If the LSP is not purely Bino but rather mixed Hig- 
gsino/Bino or Wino/Bino, the branching fraction will be 



Br. ~ 0.71 f 



(Yxg,X 
V 233 ) 



ms 



(Mixed LSP), 
(19) 

where / is the Bino fraction of the DM. 

For a Higgsino LSP the branching fraction in Eq. ( 19 ) 
will be vanishingly small since / ^ 1. However, Higgsi- 
nos can be directly produced from S decay at two-loops 
through diagrams shown in Figure [2] (middle and bottom 
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diagrams). The decay widths for the two-loop diagrams 



S^HH 



S^HH 



(hX 



1 

8^ V 471-2 
1 fhYWi 



Vt 

47r2 



mS ; 



Btt 
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^ h9i 
167r2 



( !M 



mS ; 



(20) 



where j/t ~ 1 is the top Yukawa couphng, At is the coefh- 
cient of Nf^X term in Eq. ( 10 ), and Yh = 1 is the Higgs 
hypercharge. 



Unless At is very small, the first expression in Eq. ( 20 ) 



will be dominant. Then the branching fraction for Hig- 
gsino DM will be 



Br, 



1 



167r4 



ms 



(Higgsino LSP). (21) 



Much smaller values of Br^ can be obtained in the case 
of Higgsino DM. For example, we can find Br^ ~ - i n-io _ 



10 ^ for TO5 
0.1 - 1. 



1 TeV, m 



H 



100 



10- 

500 GeV, and At 



We note that the second expression in Eq. (20) takes 



over for At ^ 0.1, which provides a lower bound on Br^^ 
regardless of how small At is. 



1^- 




H 



B H 
,AAAAA,T « 



B 



B. Higgsino Dark Matter 

In Table |lj we show the model parameters and typical 
mass scales considered in this work. In Table [TTl we give 
a summary of the various possibilities for Higgsino and 
Bino DM. 

For a purely Higgsino DM, Br^ is given by Eq. (21 1. 
Since Higgsino LSP is thermally underproduced, the cor- 
rect relic abundance can be obtained via both of the "An- 
nihilation" and "Branching" scenarios discussed in Sec- 
tion II. The important factor in determining which of the 
scenarios can work is the size of coupling At . 

To illustrate this, we first consider the case with At ^ 1. 
Then Eq. ^ resuhs in Br^ > 10"^ for > 100 

GeV. Combined with the condition for not overproduc- 



ing gravitinos from S decay, see Eq. (16 1, this leads to 
Br^ > 10~^°. Therefore the observed relic density 
cannot be obtained in the "Branching Scenario" . This 
leaves the "Annihilation Scenario" as the only possibility 
when At is large. This is shown in Table [TTl For -^t ~ Ij 
"Annihilation Scenario" is the only option even for Tj. 
at the BBN bound. The correct relic density may be 
obtained in "Annihilation Scenario" for ~ 0.5 GeV 

Next, we consider the case with At ~ 0.2. Then we see 
from Eq. (pTl) that Br^ ~ 10"^- 10"^ for ~ 100-500 



GeV. ThisTmplies that the correct relic abundance via 
the "Branching Scenario" can be obtained for 10 MeV < 
Tj. < 1 GeV. A typical example of this scenario is shown 
in Table HD 



FIG. 2: Dominant diagrams for the decay of S to Bino (top) 
and Higgsino (middle and bottom). There are additional di- 
agrams that are obtained by switching internal lines to their 
SUSY partners. 



The flexibility that one can make either of the "Anni- 
hilation" and "Branching" scenarios work by dialing At 
is also important from another point of view. "Annihila- 
tion Scenario" relies on a large annihilation rate ((Tann^^) ■ 
Therefore, tightening of the Fermi bounds on (fTannw) 
|16) . will put increasing pressure on this scenario. It is 
therefore important that one can have a viable "Branch- 
ing Scenario" , regardless of the improving constraints, by 
lowering At. 



TABLE I: Typical scales and values of model parameters in 
the scenarios considered in this work. The reheat temperature 
is Tr ~ C(l) GeV for the values presented in the Table. 



Parameter 


Value 


ms 


C»(TeV) 


mxi ~ 


0(10 TeV) 




sub- TeV 


A' 


0(1) 


At 


0(0.1 - 1) 


h 


0(10-'^) 
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TABLE II: Summary of the various scenarios considered in 
this work, for typical values of model parameters. We have 
taken ms = 1 TeV. For Higgsino DM, both "Annihilation Sce- 
nario" and "Branching Scenario" may be obtained depending 
on the value of At. For Bino DM, successful "Branching Sce- 
nario" requires light DM and low reheat, which may be ob- 
tained by suitably choosing h/mx- AH mass scales are shown 
in GeV. 



DM Mass 


Tr 


Brx 


Scenario 




100 


0.5 


6- 10^^ 


Annihilation (Aj 




H 










100 


0.5 


1 ■ 10"** 


Branching (At — 


0.2) 


B 60 3 


■ 10"^ 


4- 10"^ 


Branching (-^ ~ 


10-") 



where A„ is the coefhcient of Nu'^X term. We have cho- 
sen Alt > Ac,t, so that the loop containing u'^ dominates 
the process. We note that S decay to A'2 and quarks is 
subdominant due to the helicity suppression of fermionic 
final states, while decay to scalar N and squarks is kine- 
matically blocked, see Eq. (12 1 and the subsequent dis- 
cussion. 

A^i decay via the diagram in Fig. [3] generates baryon 
asymmetry, whose density is given by 



m = = Ys BrAT e 



(23) 



where Brjv is the branching fraction for producing A'^i 
from S decay, and e is the asymmetry parameter in Ni 
decay. Br a? is given by 



Br 



X \ 4 / X 2 



33 



■ms 



(24) 



C. Bino Dark Matter 

Bino DM is thermally overproduced. This, as pointed 
before, leaves "Branching Scenario" as the only possibil- 
ity for producing the correct DM density. In this case 
(n^/s) = Ys Br^, where Is > 2 x 10"^ for ms - 1 TeV. 
Obtaining the observed relic density then req uires that 
Br;^ < 2.5 X IQ-'^/nig. We see from Eq. ([l8| that this 
can be found for a Bino on the lighter side ~ 60 GeV 
when the reheat temperature is close to its lower value 
from BBN bound ~ 3 MeV. 

Having a viable scenario for non-thermal Bino DM is 
important. With the Fermi bounds on the annihilation 
rate improving [5^, Higgsino DM may be ruled out spe- 
cially at the lower end of the sub-TeV mass range. This 
will motivate Bino DM, for which one can find a success- 
ful non-thermal scenario based on S decay. 



BARYOGENESIS 



For \u - 0.5, we find BrAr^ - 0.01. 

Taking ~ "^1X2^ similar masses for X1.2 fermions 
and scalars. A' ~ for all flavors and colors (in ac- 

cordance with Table I), and CP violating phases of 0(1), 
the asymmetry parameter is given by 



1 

8^ 



\l X 0(10) 



(25) 



The 0(10) factor comes from summing over three flavors 
of quarks plus including diagrams where Xi are switched 
X2. Then, for A„ ~ 0.5_, we get e - 10"^ . 
It is seen from Eqs 



(|23 

TO 



24 



25 1 that we can obtain the 
0.5-1, Mat, - 500 
- 1 GeV, and mx ~ 



for Alt ^ 
10 McV 



desired value 773 ~ 10 
GeV, ms ~ O(TeV), T, 
10 - 100 TeV. 

One comment is in order at this point. N2 fermions 
produced in Ni decay themselves to three quarks. This 
must happen well before the onset of BBN in order not 
to ruin its success. The decay rate of A'2 is given by [TS] 



The entropy released in S decay dilutes any previously 
produced baryon asymmetry, see Eq. ^ , thus necessitat- 
ing mechanisms of baryogenesis at temperatures far be- 
low the electroweak symmetry breaking scale. This can 
be achieved by introducing baryon number and CP vio- 
lating operators in a suitable extension of the SM (T71IT5] . 
Late decay of S naturally provides the necessary depar- 
ture from thermal equilibrium. 

The model given in Eq. ( 10 ) can give rise to low- 
temperature baryogenesis via N decay [T7]. The decay 
of S produces TV quanta through the one-loop diagram 
shown in Fig. |3] (top diagram). The rate for S decay to 
the heavier of the two N flavors, denoted by Ni, is 



r 



8n V47r2 



V mx 



ms 



(22) 



1 



c 



126 1927r3 



TUN, 



(26) 



where B is the B— term associated with mxXX mass 
term in Eq. ( |lO| ), and C '-^ 6 is a color multiplicity factor. 
One can see that N2 decay well before 0.1 s for ttt-at^ ^ 100 
GeV, B ~ ©(TeV), and the parameter values given in 
Table I. 

Finally, we note that the "Branching Scenario" for 
non-thermal DM production (like in the case of Bino 
LSP) provides a visible sector realization of "Cladogen- 
esis" mechanism proposed in [T3] to address the baryon- 
DM coincidence puzzle. It is clear from our discussions 
that this scenario is a considerably more flexible than 
that when modulus decay is the source of non-thermal 
DM [51IT3]. 
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FIG. 3: Production of A''i from S happens at the one-loop 
level (top diagram). TVi subsequently decays and generates 
baryon asymmetry (middle and bottom diagrams). 



VI. CONCLUSION 

In large regions of SUSY parameter space, the DM 
annihilation cross-section is larger (as for Higgsino or 
Wino LSP) or smaller (as for Bino LSP) than the nom- 
inal value of 3 X 10~^^cm^ s^^ in the thermal WIMP 
scenario. Non-thermal scenarios of DM production will 
be needed in these cases to yield the correct relic abun- 
dance, while being compatible with the bounds from in- 
direct DM searches on the DM annihilation cross section 
particularly the Fermi-LAT results. 

Motivated by these considerations, in this paper, we 
have presented a non-thermal scenario that relies on the 
visible sector of a SUSY model. In this model the late 
decay of an _R-parity even scalar field S that is a SM 
singlet (but may be charged under a higher rank gauge 
group) produces DM. S is coupled to new colored fields 
X, X with a mass relation mg mx- Assuming that 
all i?-parity odd colored fields are heavier than S, it will 
dominantly decay into gluons at the one-loop level. The 
combination of a small coupling h between S and X, X, 
the mass relation mx ^ mg, and the one-loop factor 
can naturally lead to a late decay of S that yields a low 



reheat temperature Tj.. As we saw above, one can obtain 
T, - O(GeV) for ms - 1 TeV, mx ^ 50 TeV, and 
h - 10"^. 

The branching fraction for S decay to DM particles 
Br-^ depends on the nature of the LSP. If the LSP is 
Bino, then S can decay to a pair of DM particles at the 
one-loop level, where Br^^, > 10~^ for > 60 GeV. For 
a Higgsinos LSP, the same decay occurs at the two-loop 
level yielding Br^ ~ 10" — 10~^ (the exact value de- 
pending on the model parameters) for ~ 100 — 500 
GeV. As a consequence, one can obtain the observed relic 
abundance in both larger and smaller annihilation cross- 
section regions of SUSY parameter space. 

For Higgsino DM, the correct relic density can be found 
directly from S decay if Br^ is sufficiently small. The 
DM annihilation cross-section will be irrelevant in this 
case (hence "Branching Scenario" ) . We showed a bench- 
mark point for this scenario where ^ 100 GeV and 
Ti- ~ ©(GeV). For larger values of Br^, Higgsinos will 
undergo residual annihilation upon production from S 
decay (hence "Annihilation Scenario" ) . This scenario can 
yield the correct relic abundance if the annihilation rate 
is larger than 3 x 10~^^ cm^ s~^ by a factor of (Tf/Tr). 

In the case of Bino DM, where the annihilation cross- 
section is small, "Branching Scenario" is the only option. 
The correct relic density is obtained for a lighter Bino 
< 100 GeV and smaller values of reheat temperature 
T]. ^ 3 MeV. One can reduce by lowering the coupling 
h and/or raising mx- A viable scenario for non-thermal 
Bino DM is particularly important in light the Fermi- 
LAT result, which appear to prefer smaller annihilation 
cross-sections for smaller values of DM mass. 

The above scenarios are summarized in Table Ull 

The scenario presented in this paper has two main ad- 
vantages over that using moduli decay. First, gravitino 
production from S decay is naturally suppressed by the 
virtue of S belonging to the visible sector. Second, Br^ 
can be made sufficiently small by dialing the model pa- 
rameters, which is essential for a successful realization 
of the "Branching Scenario" (the only possibility in the 
case of Bino DM). To attain these virtues in scenarios of 
non-thermal DM from moduli decay, one needs to sat- 
isfy non-trivial conditions on the Kahler geometry of the 
underlying effective supergravity theory. 

Finally, our model can also successfully generate the 
baryon asymmetry of the universe. S decay also pro- 
duces the SM singlets N at the one-loop level, whose 
subsequent decay via baryon number and CP violating 
couplings to the heavy colored states X, X can yield the 
desired asymmetry t^b ~ 10~^°. Moreover, the baryon- 
DM coincidence puzzle can be addressed in the context of 
the "Branching Scenario" , since both baryon asymmetry 
and DM abundance are directly produced from S decay 
in this case. 
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Appendix A: Calculations of 5* Decay Width 

We first discuss the decay of S into two gluons. There 
is a contribution from the loop of fermions shown in 
Fig. [l] as well as a loop of sfermions. The total decay 
width is given by [19\ 



■99 



1 1 

9\/2 87r 
3 



47r2 



ms 



X tA/2(t) + t^susyW 



(Al) 



Here, mioop is the mass of the particle in the loop. The 
amplitudes A1/2 and ^susy refer to the fermion and 
sfermion loops respectively. 

With the assumed mass conditions, one obtains the 



decay width given in Eq. (13 1, which further accounts 
for the two flavors of X. 

The decay width into Binos is given by 



where 



1 

4^ 



(A4) 



where 



A/2(r) = 2(r+(r-l)/(r))- 
-4susy(t) = - (t" - /(t)) 
/(r) = arcsin^ y/r 



(A2) 



2^2 



Ana 



Yx 
2 



mxrrig x 



[Co(mx)-2C+(mx) + C+{m^)] 



and 



4"^?oop 



In the above equation, Cq and are Veltman-Passarino 
functions 120.. With the assumed mass conditions, this 



(A3) expression is reduces to Eq. (17). 
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